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Recent experimental findings and theoretical predictions suggest that nitrogen-doped CVD-grown 
graphene may give rise to electronic band gaps due to impurity distributions which favour segregation 
on a single sublattice. Here we demonstrate theoretically that such distributions lead to more 
complex behaviour in the presence of edges, where geometry determines whether electrons in the 
sample view the impurities as a gap-opening average potential or as scatterers. Zigzag edges give rise 
to the latter case, and remove the electronic bandgaps predicted in extended graphene samples. We 
predict that such behaviour will give rise to leakage near grain boundaries with a similar geometry 
or in zigzag-edged etched devices. Furthermore, we examine the formation of one-dimensional 
metallic channels at interfaces between different sublattice domains, which should be observable 
experimentally and offer intriguing waveguiding possibilities. 


I. INTRODUCTION 

The high Fermi velocity and linear electronic disper¬ 
sion in graphene appear promising for electronic devices X 
The absence of an intrinsic band gap is a potential stum¬ 
bling block for many applications. A range of possibili¬ 
ties are being investigated to redress this. Many involve 
geometric constraints in the form of, e.g ., finite-width 
nanoribbons (GNRs)^ or periodic perforations 4 An al¬ 
ternative route is the manipulation of the atomic level 
structure. The hexagonal graphene lattice is composed of 
two intersecting triangular sublattices, A and B , shown 
by hollow and filled symbols respectively in the top pan¬ 
els of Fig [TJ The equivalence of these leads to the gapless 
band structure. A sublattice dependent potential opens 
a band gap and gives mass to the charge carriers. A pos¬ 
sible implementation is to place graphene on a substrate, 
such as hexagonal boron nitride (hBN), which offers a 
potential varying on approximately the required length 
scaled. However, the potential here is quite weak and lat¬ 
tice mismatches give rise to larger scale Moire features 

Recent experiments suggest another route to breaking 
sublattice equivalence. Nitrogen-doped graphene grown 
by chemical vapour deposition (CVD) can show unusual 
distributions of substitutional N atoms. Large domains 
are found with N atoms primarily occupying a single 
sublattice^ - — This behaviour depends on growth con¬ 
ditions and theoretical works suggest possible mecha¬ 
nisms including preferential impurity positioning rela¬ 
tive to edges during growth^ and inter-impurity inter¬ 
actions in disordered ensembles j 14 v 15 Subsequent stud¬ 
ies of N-doped graphene treated by high-temperature 
annealing^, and of graphene decorated by hydrogen 
adatoms ^ suggest that asymmetric distributions may 
also arise in other scenarios. Such doping leads to dif¬ 
ferent average potentials on each sublattice and is equiv¬ 
alent to introducing an effective mass term. Extended 
graphene sheets with sublattice-asymmetric impurity dis¬ 
tributions are predicted to display electronic and trans¬ 
port band gaps, and electron-hole asymmetry in their 


conductivity 

In this work we focus on nanoribbons with sublat¬ 
tice asymmetric doping. This is motivated both by 
the possibility of etching^ and transferring^ devices 
from doped graphene sheets and by the need to under¬ 
stand the interplay between the effective mass term in¬ 
troduced by such doping and effects induced by sym¬ 
metry breaking edges. This is important since CVD- 
grown graphene contains extended edge-like defects in 
the form of grain boundaries^ - — unlike bottom-up ap¬ 
proaches which may allow synthesis of more precise 
geometries^ We are further motivated by the strong de¬ 
pendence of GNR transport on edge geometry and im¬ 
purity distribution^ - — and by sublattice dependent fea¬ 
tures in carbon nanotnbeS i 44 i 45 We consider both arm¬ 
chair (AGNR) and zigzag (ZGNR) edged ribbons, noting 
the inbuilt sublattice asymmetry of ZGNRs due to sites 
along one edge belonging to one sublattice. Similar be¬ 
haviour to bulk graphene is found for AGNRs - namely 
reliable electronic and transport band gaps consistent 
with an average mass term model. For ZGNRs, only a 
suppression of transmission is found in the expected gap 
region and it is not accompanied by a vanishing density 
of states (DOS). In particular, strong finite DOS clusters 
remain along one ZGNR edge. This is related to the po¬ 
sition dependence of simple impurity bound states near 
zigzag edges and is captured within a coherent poten¬ 
tial approximation (CPA) model. Finally, we investigate 
interfaces between different sublattice domains and pre¬ 
dict that these should give rise to robust one-dimensional 
metallic wires embedded within the gapped system, and 
which should have features detectable by scanning tun¬ 
neling microscopy (STM). 

II. MODELS 

The electronic structure of graphene is well-described 
by a nearest-neighbour tight-binding Hamiltonian with 
a hopping integral t = — 2.7eV. The use of this model 
is validated in Appendix 0 where key features from our 
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results are reproduced using a higher-order model, We 
take \t\ as the unit of energy and include substitutional 
N dopants by a change of onsite energy A = —\t\. More 
accurate parameterisations can be achieve d 22 ! 46 ! 47 but 
the qualitative behaviour described here is reasonably in¬ 
dependent of impurity species or parameterisation. We 
will discuss the change in carrier density induced by such 
dopants at the end of Section Hill below. A general band 
dispersion is given by 

e±(k) = ^{e A +e B )±l i \l (ca ~ E) 2 + 4 ~t 2 |/(k)| 2 (1) 

where ca (cb) is the potential on the A ( B ) sublattice and 
/(k) is a term arising from the sum of Bloch phases over 
neighbouring sites. For pristine graphene, ca = £b = 0.0, 
and so e±(k) = ±t |/(k)|, which is gapless near E = 0. 
Uniformly breaking the sublattice symmetry, by setting 
£a % £b, has three effects on the bandstructure: i) a 
bandcentre shift of €A + eB , ii) a direct band gap of mag¬ 
nitude \ca — £b\ a t the Dirac points and iii) the break¬ 
ing of the band linearity due to the additive constant 
(ca — €b) 2 in the square root. The quantity \ 6A ~ eB \ is 
called a mass term, and the dispersion of electrons in the 
gapped systems is no longer linear or massless. 

Transport quantities are calculated using recursive 
Green’s function (GF) techniques!!. Semi-infinite leads 
are constructed using an efficient decimation procedure!! 
and the zero-temperature conductance is given by£! G = 
, where the transmission is calculated from T(E) = 
Tr [G r r j RG G Ti y ], where Ti(E) (■i = L,R) are the level 
width matrices and G r ^ a (E) is the retarded/advanced 
Green’s function of the device region. A configurational 
average is taken for disordered systems to discern the 
overall trends. We also examine the local density of 
states (LDOS), which at site i is given by pi(Ep) = 
— ^Im [G-(Ep)] . The GFs required here involve a double 
sweep through the device region*!! 

Effective medium models are used to analyse the 
configurationally-averaged densities of states. The use of 
the two different models below allows to isolate effects 
arising from an average disorder-induced potential or 
mass-term, and the effects of scattering from individual 
impurities. Both models employ a INN tight-binding de¬ 
scription which is perfectly periodic along the ribbon di¬ 
rection. Onsite energies within the repeated unit cell are 
determined as described below. The virtual crystal ap¬ 
proximation (VCA) ignores scattering effects and simply 
takes into account the new average potential felt by elec¬ 
trons. In practise this is done by introducing a self energy 
to shift onsite energies by cA, where c is the doping con¬ 
centration and A is the shift caused by a single dopant *!! 
For sublattice dependent doping, this is generalised so 
that the self energy is sublattice dependent, Y> x = c x A x 
for x = A, B, due to c x (and/or A^) taking different val¬ 
ues on each sublattice. This new unit cell is then consid¬ 
ered part of an infinite perfectly periodic virtual crystal 
allowing us to calculate the Green’s function and thus the 


density of states. The coherent potential approximation 
(CPA) replaces this potential with a position and energy 
dependent self energy to include simple scattering effects. 
This self energy is found from the solution of the self- 
consistent equation Y> x = c x A x (l — (A x — E a ,)Geff) _1 , 
where G e ff is the Green’s function of the new effective 
medium i 51 ! 52 Periodicity of the effective medium along 
the ribbon direction can again be used to quickly calcu¬ 
late the Green’s function and density of states. It can be 
shown that the CPA includes simple scattering effects, 
but neglects higher-order scattering terms. In this way 
features appearing in the CPA, but not the VC A, arise 
due to the scattering effects beyond an average potential 
but below higher-order cluster effects such as localization. 


III. RESULTS AND DISCUSSION 

We first calculate the transmission through both GNR 
types for two disorder types - a completely random dis¬ 
tribution of impurities over all sites ( symmetric ) or a 
distribution confined to only one sublattice ( asymmet¬ 
ric ). Fig. [1] shows transmissions through a) 101-AGNR 
(width ~ 12 nm) and b) 100-ZGNR (width ~ 21 nm) 
systems. In the absence of disorder, these ribbons are 
both metallic within the nearest-neighbour tight-binding 
approximation. Results for the initially semiconducting 
100-AGNR are shown in Appendix [Aj 

The conductance of the pristine systems is shown by 
the grey shaded areas and the averaged asymmetrically 
(symmetrically) doped systems by solid red (dashed blue) 
lines. Configurational averages over 100 instances of dis¬ 
order through device regions 40 unit cells long (17 nm 
for AGNR, 10 nm for ZGNR) are shown. Impurity con¬ 
centrations are ca = 0.05, cb = 0.05 (ca = 0.1, cb = 0.0) 
for the symmetric (asymmetric) cases, where ca/b is the 
concentration on a given sublattice. Note that the asym¬ 
metric case corresponds to a random replacement of 10% 
of sublattice A carbon atoms with nitrogen atoms within 
the disordered region, for a total nitrogen concentration 
of 5% as the B sublattice is unaltered. The total concen¬ 
tration of nitrogen is thus the same for both cases. 

For AGNRs, asymmetric disorder opens a band gap 
with sharp edges on the hole side of the spectrum, in con¬ 
trast to symmetric disorder where very little transmission 
suppression is seen. The persistence of the T = 1 plateau 
in the symmetric case has been observed previously!!. In 
general, AGNRs are more sensitive to edge disorders than 
the bulk substitutional disorder considered here j 31 ! 41 The 
transport gap for asymmetric doping has a correspond¬ 
ing electronic band gap, clearly visible in the averaged 
DOS plot in c). This shows an average over the central 
800 cells of a disordered region with total length 1000 
unit cells. The appearance of this band gap is consistent 
with the results for similarly doped extended graphene 
sheet s!!. A comparison of the numerically averaged DOS 
to results from the VC A and CPA models is shown for 
the fully asymmetric case in Fig. m- Good agreement 
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FIG. 1. (Top) Schematics of a 6-AGNR and 4-ZGNR, with 
the unit cells shown by the shaded areas and the A (B) sub¬ 
lattice sites by hollow (filled) symbols. The index counts the 
dimer lines or zigzag chains across the ribbon. Remaining 
panels show results for a 101-AGNR (left) and a 100-ZGNR 
(right), a), b) show the (averaged) transmission through pris¬ 
tine systems (grey shading) and also systems with 40 unit cells 
of sublattice asymmetric (solid red lines) and sublattice sym¬ 
metric disorder (blue dashed-dotted line). c), d) show the nu¬ 
merically averaged DOS of longer systems with corresponding 
disorder profiles, e), f) show the numerically averaged DOS 
for the fully asymmetric case (black symbols) compared to 
VC A (orange) and CPA (green, dashed) model calculations. 
The concentration of N atoms for all disordered cases is 5%. 


between the VC A and numerical results is seen within 
the gap and on the electron side, while poor agreement 
is seen on the hole side. The VC A also overestimates 
the bandgap, which is somewhat smaller than the value 
c^A = 0.1 \t | given by a uniform mass term. These 
discrepancies are almost entirely corrected by the CPA, 
where excellent agreement is seen over the entire energy 
range. 

The accuracy of the VC A at gap and electron-side en¬ 
ergies suggests that the main effect of disorder here is not 
scattering, but rather an averaged potential landscape 
with a sublattice dependent mass term. The unimpor¬ 
tance of scattering effects here is also apparent in the 
transmission shown in Fig. at), where the asymmetric 
disorder only induces minor quenching of transmission 
at these energies. Conversely, the failure of the VC A 
and success of the CPA on the hole side suggest that 
scattering plays a more important role here. This is 
further evidenced by the hole-side transmission, which 


AGNR ZGNR 




FIG. 2. Transmissions for 101-AGNR (left) and 100-ZGNR 
systems with 80 unit cells of asymmetric disorder. Results 
are shown for both fully (red, solid) and partially (turquoise, 
dashed) asymmetric distributions of impurities. 


is significantly reduced relative to the pristine case and 
has its plateau features almost completely smeared out. 
This electron-hole asymmetry is consistent with results 
in graphene sheets, where reduced mobility on the hole 
side is associated with a pseudospin polarisation giving a 
higher occupation of the undoped (doped) sublattice on 
the electron (hole) side^ We have confirmed that this 
feature is also present in the AGNR case by examining 
the sublattice dependent averaged DOS. The reduced gap 
size compared to the VC A prediction is in line with a sub- 
linear gap dependence found in graphene sheets j 18 i 22 We 
have varied the concentration and find agreement with 
the Eg ~ c^( 75 scaling previously reported^ 

The right-hand side panels of Fig. [l]show that many 
of the features discussed above are radically altered for 
zigzag edged systems. Transmission suppression is ob¬ 
served in the gap region for asymmetric doping in Fig. 
CEb), but without sharply defined band gap edges. Fur¬ 
thermore, a significant DOS is noted throughout the ex¬ 
pected bandgap (Fig. [Hi). It is thus unsurprising that 
the VCA (Fig. [If) fails to capture the DOS features at 
these energies, since this model always gives a bandgap. 
However, it does capture the low-energy electron-side be¬ 
haviour quite well, including the sharp peak at E = 0. 
This peak is associated with states localised on the edge 
atoms of a ZGNR. It is doubly-degenerate in pristine rib¬ 
bons, as the states on each ribbon edge, although be¬ 
longing to opposite sublattices, are equivalent. Adding a 
uniform mass term breaks this degeneracy and the peak 
splits into two which reside at the bandgap edges, at ener¬ 
gies corresponding to the onsites of each sublattice. This 
is seen for the VCA result, but the peak at E = — 0.1 |t|, 
associated with the N-doped sublattice, is absent in the 
numerical results and only the undoped sublattice peak 
remains. The CPA once again restores the features ab¬ 
sent within the VCA, suggesting that the finite DOS in 
the expected band gap is due to scattering processes dom¬ 
inating over a gap-opening average potential. 

To verify the robustness of the gap-opening feature, 
we consider the case of less than perfect sublattice asym¬ 
metry. Fig m shows the transmissions through systems 
analogous to those in Fig. □Ji) and b), but with 75% 































4 



FIG. 3. LDOS maps of a disordered 101-AGNR (top) and 100- 
ZGNR (bottom) at E — —0.04|£|. The impurities are entirely 
on the A sublattice, corresponding to the bottom edge of the 
ZGNR, where a non-vanishing DOS is evident. 


of N atoms on sublattice A and 25% on sublattice B. 
Curves for a fully asymmetric case are shown for com¬ 
parison. For partial asymmetry, we note a clear band 
gap formation for the AGNR case, whereas transmission 
suppression without a clear band gap is still present for 
the ZGNR case. The AGNR band gap is shifted away 
from E — 0, unlike that of the perfectly asymmetric case, 
as the band centre shift and mass terms entering in Eq. 
m are no longer equal. Band gap formation at this level 
of asymmetry is promising for realising such a gap ex¬ 
perimentally, as samples with over 90% asymmetry have 
been reported^ 

To further explore the differences between armchair 
and zigzag edged geometries, we show LDOS maps for a 
single, fully-asymmetric disorder configuration of each in 
Fig EH The maps are taken at an energy in the middle of 
the expected band gap. The LDOS decays quickly as we 
move into the disordered region of an AGNR. This decay 
is also uniform across the ribbon width. For the ZGNR, 
the LDOS vanishes throughout most of the system. How¬ 
ever, large clusters of finite density remain along the 
bottom edge of the ribbon, which is associated with the 
doped sublattice. This suggests an interplay between the 
doping of a particular sublattice and the proximity of a 
zigzag edge of the same sublattice. The reproduction of 
averaged DOS features within the CPA model suggests 
that this effect can be explained in terms of single scat¬ 
tering processes, and so we now examine individual N 
dopants in a ZGNR. 

Fig 0^1 shows a few possible sites for a single N atom 
near the edge of a 50-ZGNR. The sites represented by 
red and green circles are on the edge sublattice, whereas 
the blue site is not. Fig 0}d) and c) show that impurity 
sites on the edge sublattice give rise to conductance dips 
and corresponding DOS peaks in the low energy win¬ 
dow shown here. These features, associated with anti¬ 
resonances formed by the impurity, have been studied 
previously in GNR a 3Q i 41 . Symmetry-breaking edges re¬ 
sult in a strong position dependence of the anti-resonance 
energies. Interestingly, sites near a ZGNR edge and of the 
same sublattice type can give rise to features at energies 
within the expected band gap, whereas those on the op¬ 
posite sublattice (and sites in AGNRs) result in features 
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FIG. 4. The transmission (b) and averaged DOS (c) for a 50- 
ZGNR with a single N impurity located at each of the sites 
shown by the symbol of the same colour in a). d)-f) map the 
change in LDOS near the three possible impurity locations, 
taken at E — —0.05|t| (d) and E — — 0.02|t| (e and f). 


at energies far outside this window. In Fig 0] d)-f) the 
change in LDOS near three of these sites is mapped. For 
d) and f), corresponding to sites on the edge sublattice, 
we choose the DOS peak energy and note a significant 
triangular region of increased DOS near the impurity lo¬ 
cations at the bottom edge. For the opposite-sublattice 
impurity site in e), we choose the same energy as d), and 
note that no such feature is visible and the DOS barely 
differs from that of a pristine ribbon. Consequently elec¬ 
trons in this energy range are scattered by impurities 
located on the same sublattice as the edge, and not by 
those on the opposite sublattice. Returning to asymmet¬ 
rically disordered ZGNRs, we can understand the finite 
DOS in the expected bandgap (Fig [l]i) as the average 
of many single impurity peaks at different energies and 
corresponding to A-sublattice impurities at different lo¬ 
cations near the bottom edge. Away from this edge, the 
density of states vanishes as shown in Fig. [3j because 
the net effect of the doping here is an average mass term 
and not scattering from impurity states. This is con¬ 
firmed further by examining the position dependence of 
the CPA self-energy, FA, which in AGNRs takes a real 
and quite uniform value slightly smaller than c^A. This 
is also true across much of a ZGNR, except near the edge 
associated with the doped sublattice, where FA becomes 
complex and its real part varies drastically from c^A. 
The VC A is unable to explain behaviour near this edge, 
as the net effect of the doping is no longer an effective 
mass term. Increasing the device length will lead to a 
transport gap as we enter the localization regime. How¬ 
ever, this gap is unrelated to the effective mass term or 
a DOS gap, and is similar to the behaviour observed for 
ZGNRs with symmetric doping. 

The breakdown of the band gap in asymmetrically 
doped graphene near a zigzag edge may have interest¬ 
ing consequences beyond GNR devices. Grain bound¬ 
aries can have geometries similar to zigzag edges and 
break the lattice symmetry in the same manner—. Thus 
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FIG. 5. LDOS maps of asymmetrically doped 200-AGNRs 
with a) sudden or b) gradual sublattice interfaces running 
along the centre of the ribbon, taken at E = — 0.04|£|. c) 
shows the transmissions for these systems compared to one 
with a single domain. 


CVD-grown systems may experience leakage near these 
boundaries. Another relevant interface is that between 
neighbouring regions with doping on opposite sublattices. 
These have been mapped experimentall}^, and in Fig. [5] 
we consider an AGNR with a sublattice interface running 
parallel to the edge so that only the A (B) sublattice is 
doped in the bottom (top) of the device. Near the bound¬ 
ary the average mass term switches sign, closing the band 
gap and resulting in states confined near the interface^. 
This is confirmed in the LDOS maps in Fig. [5j shown 
for systems with both a ) sudden and b ) gradual interfaces 
where the impurity concentration changes linearly from 
one sublattice to the other over 4 or 20 atoms respec¬ 
tively. In both cases we note a large, finite DOS running 
along the interface and decaying away from it. Further¬ 
more this regions acts as a propagating channel as is clear 
from panel c), where a finite transmission is noted across 
the band gap region of a single domain device. 

Electron doping by nitrogen impurities shifts the Fermi 
energy, Ep, relative to any gap. Accessing the gap re¬ 
gion experimentally will involve the application of a gate 
voltage. While accurate electron-counting can be per¬ 
formed within DFT calculations^ for single impurities 
or small disordered regions, this is not feasible for the 


system sizes considered here or in experiment. Nonethe¬ 
less, the charge density fluctuation can be approximated 
from Sn ~ e d( c a+c b )pc ^ w h ere e d ~ 0.4 i s the average 
doping efficiency of nitrogen in GNRs^i and pc is the 
density of lattice sites in graphene. For ca = 0.1, we find 
Sn rsj 7.6 x 10 13 cm -2 , just inside the range of the most ad¬ 
vanced gating methods^ ca — 0.02 gives a more realistic 
Sn ~ 1.5 x 10 13 cm -2 , while yielding Eg ~ 50—200 meV. 
Gaps from lower concentrations, whilst too small for ap¬ 
plications, still allow experimental verification of our re¬ 
sults. It is also possible to shift Ep nearer the gap by 
codoping with a symmetrically distributed p-dopant, at 
the cost of reducing transmission outside the band gap. 


IV. CONCLUSIONS 


Our results highlight the importance of edge geometry 
in doped graphene nanoribbons. The band gap predicted 
for sublattice asymmetrically doped graphene is sensitive 
to the presence of zigzag edges, where a gap-opening av¬ 
erage potential is no longer the dominant effect of disor¬ 
der. Instead impurity bound states within the expected 
band gap, associated with the edge sublattice, lead to a 
finite DOS and propagation, albeit scattered, along the 
edge. A band gap opening, similar to that in graphene 
sheets, is observed for armchair edges. The sensitivity of 
gap opening to edge geometry is relevant beyond ribbon 
devices. The majority of samples with sublattice asym¬ 
metric disorder are grown by CVD, which gives rise to 
edge-like defects in the form of grain boundaries. Since 
these can have zigzag-edge like symmetries, we expect 
similar leakage near grain boundaries in asymmetrically 
doped poly crystalline graphene sheets. This may make 
it difficult to verify experimentally the band gaps pre¬ 
dicted for such systems. Finally, we show the formation 
of one-dimensional metallic wire behaviour at the inter¬ 
face between two regions with doping on opposite sub¬ 
lattices. Such interfaces are present in experimental sys¬ 
tems, and the features we predict should be observable to 
STM measurements. These channels present waveguid- 
ing possibilities as, away from defects or edges, leakage 
is prevented by the gapped region surrounding them. 
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FIG. 6. Transmissions for pristine (top) and asymmetrically 
disordered (bottom) ribbons using both INN (solid, black 
curves) and 3NN (red, dashed curves) models. The INN re¬ 
sults for the 101-AGNR and 100-ZGNR are reproduced from 
the main text, whereas the 100-AGNR case represents an ini¬ 
tially semiconducting ribbon within the INN model. 


Appendix A: Comparison of 1st and 3rd 
nearest-neighbour tight-binding results 

To check the validity of the first nearest-neighbour 
tight-binding approximation (INN) for our systems, we 
compare the transmissions of pristine and asymmetrically 


disordered nanoribbons calculated with both this model, 
and with a more complete third nearest-neighbour (3NN) 
description of graphene. We also consider a 100-AGNR 
which is semiconducting within a INN description in the 
absence of dopants. The INN results are based on the 
system in Figs up and b), where a constant value of 
t = — 2.7eV is used throughout the system to describe 
the hopping parameter between nearest neighbour sites. 
The 3NN results are calculated using the same relative 
2nd and 3rd neighbour hoppings for pristine graphene as 
in Ref. [22|. For both models, we use a simple onsite shift 
of A = —\t\ to represent an impurity. Larger values of 
A, suggested elsewhere in the literatur o 22 i 46 for nitrogen, 
would enhance the features discussed in this work due 
to the scaling of the effective mass term with scatterer 
strength. 

For AGNRs we note that the INN model captures all 
the main features, with the exception of the previously re¬ 
ported small band gap for pristine 101-AGNRs. We also 
note the band gap opening induced by asymmetric dis¬ 
order occurs regardless of the metallic or semiconducting 
nature of the corresponding pristine ribbon. The higher 
transmission values for pristine ZGNRs at low electron- 
side energies are due to the zero-energy peak no longer 
remaining disperionless within the 3NN model. This has 
been reported previously in the literature^ We note that 
the key result discussed in our paper, namely the band 
gap opening or transmission suppression at low hole-side 
energies in asymmetrically doped systems, are perfectly 
described by the INN model. 
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